CMB anisotropies from primordial inhomogeneous magnetic fields 
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Primordial inhomogeneous magnetic fields of the right strength can leave a signature on the 
CMB temperature anisotropy and polarization. Potentially observable contributions to polarization 
B-modes are generated by vorticity and gravitational waves sourced by the magnetic anisotropic 
stress. We compute the corresponding CMB transfer functions in detail including the effect of 
neutrinos. The shear rapidly causes the neutrino anisotropic stress to cancel the stress from the 
magnetic field, suppressing the production of gravitational waves and vorticity on super-horizon 
scales after neutrino decoupling. A significant large scale signal from tensor modes can only be 
produced before neutrino decoupling, and the actual amplitude is somewhat uncertain. Plausible 
values suggest primordial nearly scale invariant fields with B\ ^ 10~^''G today may be observable 
from their large scale tensor anisotropy. They can be distinguished from primordial gravitational 
waves by their non-Caussianity. Vector mode vorticity sources B-mode power on much smaller scales 
with a power spectrum somewhat similar to that expected from weak lensing, suggesting amplitudes 
Bx ^ 10~^G may be observable on small scales for a spectral index n ~ —2.9. In the appendix we 
review the covariant equations for computing the vector and tensor CMB power spectra that we 
implement numerically. 



I. INTRODUCTION 



Magnetic fields are ubiquitous in the universe, with 
~ 10~^G coherent fields observed on galactic and cluster 
scales. However their origin is not well understood (for 
a review see Ref. [l|). Tiny seed fields <■ 10^^°G may 
have been amplified by a dynamo mechanism to give the 
much larger fields we now see, though to what extent this 
process can work in practice is not yet clear 2]. Alter- 
natively initial fields of strength ^ 10~^G can give rise 
to galactic fields of the observed values without a func- 
tioning dynamo mechanism. Such fields have potentially 
interesting observational signatures on the CMB, and if 
present would provide powerful constraints on models of 
the early universe. The absence of such signatures may 
also serve as a consistency check on models of galaxy evo- 
lution that would be observationally incompatible with 
initial fields this large. 

A primordial field of ~ lO^^G can leave a signature 
in the B-mode (curl-like) CMB polarization. Since the 
scalar (density) perturbations do not produce B-modes at 
linear order, the B-modes are a much cleaner signal of ad- 
ditional physics than very small fractional changes to the 
temperature or E-mode polarization. However B-modes 
are produced at second order through weak lensing 3, j|| , 
and are also generated by primordial gravitational waves 
(tensor modes) M ■ Other possible sources include topo- 
logical defects 0, H- The focus of many future CMB 
observations will be on observing the B-modes, so it is 
useful to assess in detail the various possible components 
and how they can be distinguished from each other. 

Primordial fields with a blue spectrum compatible with 
nucleosynthesis are far too weak on cosmological scales to 
leave an interesting signature Ql- In this paper we con- 
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sider in detail the CMB signal expected from ^^ lO^^G 
primordial fields with a nearly scale invariant spectrum. 
Such fields are not well motivated by current theoreti- 
cal models, which mostly give much smaller amplitudes 
or a much bluer spectrum [ll M> HU UM ■ Observation of 
primordial fields at this level would therefore be a pow- 
erful way to rule out many models. However Ref. p^ 
present one model in which observably interesting CMB 
signatures may be produced. 

Previous semi-analytical work has investigated the 
CMB signal from both tensor [ij, [lj| and vector 
modes [13, u^ sourced by magnetic fields. Here we give 
a more detailed numerical analysis of the full linearized 
equations. We include the effect of neutrinos as they 
change the way that magnetic fields source gravitational 
waves and vorticity on super-horizon scales. Our final 
CMB power spectra include the contribution to the B- 
mode signal from both the tensor and vector modes. We 
do not consider helical modes [la, [13 which can be de- 
tected via their parity-violatin g co rrelations, nor the ef- 
fects of Faraday rotation [l^, [l3 (which can be iden- 
tified by the frequency dependence). We also assume 
that reionization is relatively sharp and unmodified by 
energy injection into the IGM from decay of the small 
scale magnetic field [23 . For a discussion of constraints 
on homogeneous magnetic fields see e.g. Ref. [2l| and 
references therein. 



II. COVARIANT EQUATIONS 

We consider linear perturbations in a flat FRW uni- 
verse evolving according to general relativity with a cos- 
mological constant and cold dark matter, and approx- 
imate the neutrinos as massless. Perturbations can be 
described covariantly in terms of a 3-f 1 decomposition 
with respect to some choice of observer velocity Ua (we 
use natural units, and the signature where Uau" = 1), 



following Refs. |2a, |2a IM |25| • The stress-energy tensor 
can be decomposed with respect to Ua as 



Tab = PUaUb - phab + 2U(^aqb) + T^ab 



(1) 



where p is the energy density, p is the pressure, qa is the 



heat flux and and 



T^ab 



T(ab) is the anisotropic stress. 



Angle brackets around indices denote the projected (or- 
thogonal to Ua) symmetric trace-tree part (PSTF). The 
tensor 



Qij "^iUj^ 



(2) 



where Qij is the metric tensor, projects into the instan- 
taneous rest space orthogonal to Ua- It defines a spatial 
derivative D^ = ha^^b orthogonal to Ua where Va is 
the covariant derivative. Spatial derivatives can be used 
to quantify perturbations to background quantities, for 
example the pressure perturbation can be described co- 
variantly in terms of Y>aP- 

Conservation of total stress-energy V°'Tab — implies 
an evolution equation for the total heat flux qa 



9a + o®9a + {P+P)Aa 



DaP 



D^TT, 



ah 







(3) 



where qa = u^^bqa, © = V'^Ua is three times the Hubble 
expansion, and Aa = Ub^^Ua is the acceleration. The 
evolution equation for the heat flux g^ = (p' + p'^)v\ of 
each matter component present is of the form 



so the magnetic field simply redshifts as l/S"^ where S 
is the scale factor, and hence tt^j, oc p^ (x 1/5"^. More 
general equations can be found in e.g. Ref. |25| . 

The Poynting vector heat flux is zero {q^ = 0) in all 
linear frames because we have set Ea = 0. Since Aa is 
first order, on linearizing we have the Lorentz force L^ 
given by the evolution equation Q 



^b_-B 



D^7r-,^DaP"=L-. 



(9) 



This is consistent with the usual curlB x B expression. 
The opposite force acts on the baryons to ensure total 
momentum conservation, which with the Thomson scat- 
tering terms |24| gives the baryon velocity evolution equa- 
tion: 



Va + -^QVa +Aa ^ 

3 p" 



p 



,7 



UeCTT 1 -^Va - la 



P'' 



(10) 



where rie is the electron number density, ax is the Thom- 
son scattering cross-section, p'^ is the photon energy den- 
sity, and we neglect baryon pressure terms of the form 
p'^ I (}' <C 1. Thus magnetic fields source baryon vor- 
ticity, as well as providing extra density and pressure 
perturbations, and anisotropic stresses. 

We define the vorticity vector fia = curlua where for 
a general tensor 



9a + 369! -f (p' 



-p')Aa-Y)aP' 



^'<b = LI 



(4) 



where L^ is an interaction force term. Conservation of 
total stress energy qa = X^i Qa implies that ^- L\ = 0. 
For magnetic fields the components of the stress-energy 
tensor are given by^ 



P 



T^ab 



'ip' 



i(^^ 



B' 



-E{aEb) — Bl^aBb) 

-{ExB)a 



(5) 

(6) 

(7) 



where Ea and Ba are the electric and magnetic field pro- 
jected vectors. We take B^ and E'^ to be first order, and 
to this order Ea and Ba are frame invariant. For most 
of its evolution the universe is a good conductor so we 
may take Ea — in all linear frames: the magnetic fields 
are frozen in, and in this approximation almost all the 
complications of MHD disappear. The linearized Bianchi 
identity for the electromagnetic field tensor implies 



Ba + -QBa = 



(8) 



curl Xa^...ai 



7^,,a^a^u'D-X^ 



a2...ai) 



(11) 



and round brackets denote symmetrization. It is trans- 
verse D°f2a = 0. Remaining quantities we shall need are^ 
the 'electric' Eab and 'magnetic' Hab parts of the Weyl 
tensor Cabcd 

Eab = CacbdU^u'^ Hab = -jV acdf G be"^ u"" U^ (12) 

(which are frame invariant) and the shear Uab = D/a^h). 
The Einstein equation and the Bianchi identity give the 
constraint equations 






1 



■ curl rih 



-Dfc8 — Kf/b — 



k( ^gh+gObp+iov^b) =0 



1 



n[{p + p)rib+ curl qb] = 



Ha 



curl fJn 



1, 



-D/„fi, 



{a^l'b) 



(13) 



^ Note that unlike many other authors we use natural rather than 
Gaussian units. Due to the signature choice —E^ > 0. 



^ From here on we do not use E for the electric field; E^b has 
nothing to do with electromagnetism, and is merely called the 
analogous 'electric' part of the Weyl tensor by analogy. The 
'electric' part of the polarization distribution is written as £^j . 



and the evolution equations 



III. TENSORS 



^a + o©^o = CUrMa 



CTah + 7:^'^ab 



Eab + 0-Eab 
Hab + ©-ffab 



— Eab — -zKT^ab 
curl i/„(, + 



TTab - (/O + p)(Jab + I^T^ab 



curl Eab - -^CUrlTTafc. 



(14) 



Here k = 8ttG. 

The distribution functions for the various species can 
be expanded into multipole moments. For example the 
photon multipole tensors Ia, = I{ai...ai) are defined as 
moments of the distribution of the photon intensity /(e) 
per solid angle as [2g 



Ia, 



dVteI{e) e(A^), 



(15) 



where the direction vector e^ is normalized to e^Ca 



-1 



and 



HAi) 



3(ai 



°o.i) 



are irreducible PSTF tensors. 



The ei^Ai) ^^^ orthogonal: 



— / dVL, 

An 



HAi) 



.(B„)_,, t&l.i"^ 



{21 + 1) 



I (ai 



ai) 



(16) 



The Iai multipole tensors have 21 + 1 degrees of freedom, 
lab = T^Zb ^^ ^^^ anisotropic stress, la — q2 is the heat 
flux and I = p-y. The temperature anisotropy can then 
be expanded as 



AT(e) 



E 



(2^ + 1)! lA,e' 
A{-2)\l\Y p. 



EE 

/ m— — / 



(^ImXhni^) 



(17) 
where the latter expansion in terms of spherical harmon- 
ics Yim is the non-covariant version of the expansion in 
Iai ■ The CMB power spectrum is defined in terms of the 
variance of the spherical harmonic components a/,„ by 



Ci = 



\air. 



IT {21)1 {IaJ^') 
4{-2)^{U) 



P'. 



(18) 



Analogous results for the polarization are given in 
Ref. 27], where £ai is a gradient-like multipole of the 
polarization tensor and Bai is a curl-like multipole. 

The covariant equations can be expanded in terms of 
scalar, vector and tensor harmonics. The details and 
definitions are given in the appendix. In the follow- 
ing sections we analyse in detail the tensor and vector 
equations, where each quantity is a component of a har- 
monic expansion, and k labels are suppressed. Scalar 
modes can source temperature and E-polarization CMB 
signals, however since we are mostly interested in the B- 
polarization signal, which is not sourced by scalar modes, 
we do not discuss scalar modes here. A partial analysis 
of scalar modes is given in Ref. I28j . 



Expanding in ttt, = 2 tensor harmonics (|A12|1 (and sup- 
pressing m indices), the constraint equations (|13l) imply 
that the Weyl tensor variable H is related to the shear 
hy H = a. The evolution equations (|14|l then give 

k^{E' + HE)-k^a+'^S^{p + p)k(7 = '^S^{n' + Hn) 



a' + 21-1(7 + kE 



2k 



(19) 



where primes denote derivatives with respect to confor- 
mal time 77 and Ti = 59/3 is the conformal Hubble pa- 
rameter. The Weyl tensor variable E and the shear a 
define the new variable 



Ht = -2E - 



k 



(20) 



to correspond to the metric perturbation variable of non- 
covariant approaches. It satisfies Hip = —ka, and the 
above equations combine to give the well known evolution 
equation 



Urp 



2HHL 



k^Hr = nS'^n. 



(21) 



Magnetic fields provide a component of the anisotropic 
stress n and hence source gravitational waves, and we 
quantify the magnetic field source by the dimcnsionless 
ratio -Bo = H^/p-y. The covariant tensor equations are 
discussed in more detail in Ref. |2g. 

Equations for the evolution of the tensor multipoles are 
obtained from the appendix {m = 2 in Eqs. IjAlSp and 
(|A14p ) . We use a series expansion in conformal time 77 to 
identify the regular primordial modes in the early radi- 
ation dominated era. Defining lu = ftm'Ho/V^Rj where 
riji = Q^ + fl^, and Tip and Qi arc the Hubble parameter 
and densities (in units of the critical density) today, the 
Fricdmann equation gives 



S 



^ 'm "0 



UJ^ 



uj'q+-uj^T]'^ + 0{t]^) 



(22) 



Defining the ratios R^ = il.^/flR, R^ = fl~f/flfj and keep- 
ing lowest order terms the regular solution (with zero 
initial anisotropics for I > 2) is 



Hr, 



Hr 



(0) 



5 (A:r/)2 



(0), 



4i?^ + 15 
R-fBo 



Rv 



2 4i?^ + 15, 

15 R-fB^kr] 
' 14 4i?^ + 15 
, 15 (fcr;)2 



15 R^Bo{kf])^ 
28 4i?^ + 15 



(23) 



i{krffH^ 
UiR^ + ldJ '" 3 4i?^ + 15 



r(o) 



where H}p is the initial value (after neutrino decoupling) 

and TT^ = H^/p^. The Bq 7^ mode (with H^'^ = 0) 
has compensating anisotropic stresses: the sum of the 



kS'^h 



45 R^k^Bo 
14 4i?^ + 15 



(45 - 2R^)Lor] 
2R^ + 15 



Oiv^ 



magnetic and neutrino terms gives the total source term Since the transverse traceless part of the metric tensor 

hij — ^^ _j_ 2HtQ^ij this corresponds to a primordial 
power spectrum for h of 

Ph^[m^Hr^':imn)fPB,- (29) 

Power spectra are defined in Eq. HA26|I and R-y ^ 0.6. 

Thus the tensor covariance part of the magnetic field 
signal is very similar to that expected from primordial 
gravitational waves, and can be computed trivially by 
using the above tensor power spectrum in the numerical 
codes CAMB "2^ or CMBFAST [s^. However unlike 
primordial gravitational waves the spectral index of Ph 
here is expected to be at least slightly blue, and the signal 
should be non-Gaussian because Bq is quadratic in the 
magnetic field. Allowing for subtracting the B-mode lens- 
ing signal, levels of Ph ~ 10~^^ may be observable'' [3ll |. 
This corresponds to Bx ~ 10-i°G for rjir./v'i ~ 10"® and 
a close to scale invariant spectrum. To distinguish this 
from primordial gravitational waves from inflation one 
would need to detect the non-Gaussianity or small scale 
power from the vector modes. 

The compensation mechanism in principle works with 
any coUisionless relativistic fluid, even when it only 
makes up a fraction R^, ^ of the energy density [R^ 
can be interpreted as any coUisionless component). How- 
ever if there is coUisionless relativistic matter only from a 
time 77* after magnetic field production, partial compen- 
sation will be effective at a time 77 ~ 77*6^"'/^" . For small 
fractions Ry this is a very large time, so the mechanism is 
inefficient for components that are only a small fraction of 
the density allowed by the nucleosynthesis bound. Neu- 
trinos themselves can only suppress gravitational wave 
production after neutrino decoupling at z ~ 10^. 



(24) 

rather than the oc 5^p^ ex 1/77^ result expected without 
coUisionless radiation. For fc <C 7i there is therefore no 
sourcing of gravitational waves during radiation domina- 
tion, so Ht oa {kriY' if it was zero initially. CoUision- 
less fluids suppress generation of gravitational waves on 
super-horizon scales. 

Before neutrino decoupling there is no neutrino 
anisotropic stress so the magnetic field source is not com- 
pensated. Taking rji^ as the magnetic field production 
time (at which we take Ht = cr = 0), the solution for 
^77 ^ 1 is approximately"^ 



(25) 
(26) 



Ht « 3i?^Bo ln(77/7/i„) + ^ - 1 



■n 



SR'jBq 
krj 



1 



V 



After neutrino decoupling this mode must convert into a 
combination of the above regular modes, which can then 
be used to compute the observable signature. As the 
neutrino coupling is switched off the neutrino anisotropic 
stress becomes important, and with no scattering evolves 
as 






15 



rkc 



(27) 



Since the octopole J3 and tt^ will be zero before neutrino 
decoupling, for modes well outside the horizon we have 
ttI ~ —Bq/t] just after decoupling (we assume 77 ^ rjin), 
and hence the neutrino anisotropic stress grows logarith- 
mically TTy ^ ^Bq ln{rj/rj^). It therefore reaches the con- 
stant value TTi, ^ —Bq of the regular solution in about 
one e-folding. At this point Ht ceases to grow logarith- 
mically because the magnetic anisotropic stress source is 
now cancelled by the neutrinos, and Ht gives the ampli- 
tude of the usual regular solution Hj. . 

Thus after neutrino decoupling we expect a combina- 
tion of the usual passive primordial tensor mode and the 
regular compensated sourced mode. As we show explic- 
itly below, the compensated mode can be neglected com- 
pared to the small scale vector mode contribution. The 
passive tensor mode has 



IV. VECTORS 

Unlike tensors modes, vectors are not in general frame 
invariant. It is therefore convenient to choose the frame 
Ua to simplify the analysis. At linear order one can al- 
ways write Ua — u^ + Va , where w^ is hypersurface or- 
thogonal and Va is first order, so cwclua — curlwa. For 
a zero order scalar quantity X it follows that 



D„X = B^X 



,X. 



(30) 



H 



(0) 



3R^Bo ln(77^/77i„) 



(28) 



where rj'^ is the time of neutrino decoupling (assum- 
ing magnetic field generation is during radiation dom- 
ination at 77 ~ 77i„ and before neutrino decoupling). 



For vector modes (D^X)*^^) = 0, and it is convenient 
to choose the frame Ua to be hypersurface orthogonal so 
that curlua — and hence (DaXY^^ = 0, where the bar 
denotes evaluation in the zero vorticity frame. From its 



It appears that a will be large when fcr; <^ 1, however the con- 
tribution to (a"a6f"')/(Kp7) and similar terms remain small, so 
this docs not violate the linearity assumption. 



* Observable in the sense that in the null hypothesis that there 
are no non-lensing B-modes, the residual noise level would be 
consistent with this level. Actual subtraction of CMB lensing 
in the presence of B-modes from magnetic field sourced vector 
modes may be extremely difficult, but one should still be able to 
detect violations of the null hypothesis that there are none. 



propagation equation, vanishing of the vorticity also im- 
plies that Aa = 0, so the zero vorticity frame coincides 
with the synchronous gauge. The CDM velocity is also 
zero in this frame modulo a mode that decays as 1/ S 
where S is the scale factor. 

Expanding in the m — 1 vector harmonics ljA12|l the 
equations for the harmonic coefhcients in the zero vortic- 
ity frame (fi = 0) reduce to 



k{a' + 2na) = -nS^U 



H = -a 

2 



2KS^q^Pa. 



(31) 



The combination v + a (the Newtonian gauge velocity) 
is frame invariant, as are a = a + ft and v = v — fl. By 
choosing to consider the zero vorticity frame we have sim- 
ply expedited the derivation of the above frame invariant 
equations. Other papers use the Newtonian gauge |3a |. 
in which a is the vorticity. The equations are consistent. 
The baryon vorticity evolution equation H10() becomes 



v' + nv 



Pb 



SUeCTT I -rw 



h 



1 



-kBn 



(32) 



where /i = 4w-y/3 and as before _Bo = I\-b/ P-y- 

At early times the baryons and photons are tightly 
coupled, the opacity t^^ = Sn^aT is large. This means 
w^ « w, and we can do an expansion in Tc that is valid 
for e = max(A:Tc,7iTc) <SC 1. To lowest order the baryon 
velocity evolves as 



R 



l + R 



Hv 



3kBo 



+ 0(rc) (33) 



where R = 3pf,/Ap-y. 
ticity is 



The solution with zero initial vor- 



3 Bokrj 
"Sl + R' 



(34) 



so the magnetic field sources a growing baryon vortic- 
ity. At matter domination the vorticity starts to redshift 
away, however by recombination this will only be an or- 
der unity effect. On smaller scales where fcT^ = 0(1) be- 
fore decoupling the perturbations are damped by photon 
diffusion, giving a characteristic fall off in perturbation 
power on small scales. 

To identify the primordial regular mode we perform a 
series expansion as we did for the tensor case. Assum- 
ing no primordial radiation vorticity (the regular vector 
mode [33) the result is 
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(35) 
(36) 
(37) 
(38) 



As in the tensor case the anisotropic stresses compensate 
on super-horizon scales, so they source negligible shear 
a on these scales. However, unlike the tensor case, any 
non-zero a present on super-horizon scales at neutrino 
decoupling decays away rapidly and has no observable 
effect, so the evolution prior to neutrino decoupling is 
irrelevant. The observable signature comes from the vor- 
ticity sourced on sub-horizon scales by the magnetic ani- 
siotropic stress on its own. 

In the approximation that recombination is sharp at 
7] — rj* , the photon multipoles are given approximately 
from the integral solution ljA16|l by 



him) 



(w-hcr)*; -I- 



4 dx 



dr^a'^i, (39) 



where '^i = lji{x)/Xj ji{^) is a spherical Bessel function 
and X = ^(^0 ~ v)- Thus the Newtonian gauge vorticity 
a has last scattering and integrated Saches- Wolfe (ISW) 
contributions to the temperature anisotropy. In the ab- 
sence of neutrino anisotropic stress a ~ — 3i?^-Bo/(fc?y) 
during radiation domination, as in the tensor case. On 
super-horizon scales this is large, and would give a large 
scale contribution from a orders of magnitude larger 
than the small super-horizon Doppler contribution from 
V. Previous work [Ij, |3J| has neglected the contribu- 
tions from a around last scattering, an approximation 
that is good on small scales. However on super-horizon 
scales when a is large, the diverging a contributions are 
far from negligible. Previous power spectra are approxi- 
mately the right shape on large scales, but for the wrong 
reason: the large scale anisotropics are small because of 
neutrino compensation suppressing the source term for ct, 
not because they are insensitive to a. Similar comments 
apply to the polarization power spectra. 

There is one caveat to the above. The decay of a from 
neutrino decoupling to last scattering amounts to a decay 
factor of about (z^/z'J)'^ ^ 10^^. However on arbitrarily 
large scales the l/{krf^) evolution of a before neutrino 
decoupling can be larger than this, so on the very largest 
scales there can still be a contribution to a at last scat- 
tering. In the asymptotic limit the dipole /i appears 
singular, though the quadrupole I2 is finite. The con- 
straint n > —3 ensures that the total power from large 
scale modes is not singular if the individual modes are 
not. Here we neglect this effect, effectively assuming the 
power spectrum cuts off on sufficiently large scales that 
are otherwise unobservable. It is unclear whether there 
is a serious infrared problem with a nearly scale invariant 
spectrum or not. We simply compute the CMB transfer 
function from a given anisotropy stress Bq, and defer the 
issue of what the spectrum actually is, how it could be 
generated, and its actual asymptotic behaviour. 



V. NUMERICAL RESULTS 

In this paper the focus is on calculation of accurate 
CMB transfer functions from a given initial distribution. 





FIG. 1: Typical CMB temperature TT (top solid), po- 
larization EE (bottom solid), BB (dashed thick) and cross- 
correlation TE (dotted; absolute value) power spectra from 
vector modes with Bx = 3 x 10~®G, n = —2.9. The thin 
dashed line shows the scalar adiabatic mode TT spectrum 
(without magnetic fields). The increase in power at £ < 10 is 
due to reionization at redshift z ~ 13. 



As a convenient ansatz for computing sample Ci power 
spectra we assume a Gaussian primordial magnetic field 
distribution, with power spectrum Pb oc k"~^ (the defi- 
nition of n is conventional). 

One can define a smoothed magnetic field B\ using a 
Gaussian smoothing of width A (we choose A = IMpc) 
and express the power spectrum in terms of B\ as in 
Refs. ji3;llJ|- ^^ harmonic space the anisotropic stress is 
given by a convolution of the underlying magnetic fields, 
so the power spectrum for the anisotropic stress at a given 
k feels the power from across the Pb spectrum. For vec- 
tors and tensors the resultant power spectrum is given 
approximately by jl4| 



(27r 



|"+3B2 



(2n + 3) 

, \ 2n+3 



2r(^)p, 



3 



2n+6 



(40) 



for — 3 < n. The scale kjj comes from a small scale damp- 
ing cut-off Tj?, iA\ , and does not affect the power spec- 
trum significantly for nearly scale invariant power spectra 
with n ~ — 3. The spectrum is singular at n = — 3, which 
comes from the singular build up of super-horizon power 
for a scale invariant B spectrum with no large scale cut- 
off. Since Bq = H/p-y is quadratic in B, the spectrum 
of Bq will be non-Gaussian, so the power spectrum only 
contains a subset of the available information, though it 
is useful to assess the detectability amplitude. 

As a sample example we take B\ = 3 x lO^^G and 



FIG. 2: Typical CMB B-mode polarization power spec- 
tra spectra for vector modes (thick dashed), tensor modes 
(dotted) and total (solid) for Bx = 3 x lO'^'G, n = -2.9, 
Vt/Vin ~ 10®. The bottom dotted line is from the negligible 
tensor modes sourced after neutrino decoupling (the compen- 
sated mode). The top dotted line is from tensors sourced af- 
ter magnetic field generation until neutrino decoupling (and 
should be regarded as an estimate correct to a few orders of 
magnitude). The thin dashed lines show the B-mode spec- 
trum from weak lensing (peaking at ^ ~ 1000), and scale 
invariant primordial tensors with initial power ratio ~ 10~^ 
(peaking at ^ ~ 100). The magnetic field spectra scale as B^- 



n = —2.9 (as in Ref. J5j), which implies 

0.2 



P, 







1.16 X 10"i3 ( T^ 



(41) 



Since data will only constrain Pbo directly, we take this 
equation to be exact for our numerical results so they 
may easily be related to different power spectra for Pb^ 
(which may or may not come from the assumed power 
law spectrum of Ba fluctuations). The power spectrum 
Pbo scales as B'l as do the CMB power spectra, so large 
changes in overall amplitude can be obtained from rela- 
tively small changes in the primordial field: the value of 
B\ has to be quite finely tuned to give a CMB signature 
that is neither totally dominant nor totally negligible. 

Numerical results from vector modes with B\ = 3 x 
10~^G are shown in Fig.^ in comparison with the spec- 
tra expected from primordial curvature perturbations 
and possible primordial gravitational waves. For this B\ 
the effect on the temperature power spectrum is negligi- 
ble; only ifBA^SxlO^^G could there be a significant 
contribution to the power at I > 2000, perhaps contribut- 
ing some of the power observed on these scales |33 . 

The contributions to the most easily distinguished B- 
modes are shown in Fig. [3 including the tensor contribu- 
tion. It is clear that the compensated tensor mode has a 



negligible observational signature and can be neglected. 
The exact amplitude of the large scale tensor signal from 
gravitational waves sourced before neutrino decoupling is 
uncertain because we do not know the time (or mecha- 
nism) of field generation, nor have we modelled neutrino 
decoupling in detail. 

Our vector mode results are in broad agreement with 
the semi-analytical results of [15] . The main qualitative 
difference is that our TE cross-correlation changes sign in 
the damping region. The quantitative results differ some- 
what across the spectrum due to our more detailed full 
analysis of the damping, recombination, inclusion of neu- 
trinos and modelling of reionization (we have also used 
a slightly different primordial power spectrum). The re- 
sults in Ref. [1J| for the large scale vector and tensor spec- 
tra are too large by a factor^ of {2tt)^ (giving constraints 
on Bx too small by a factor of (27r)'^/'* ^ 4). There was 
another normalization error in 36] but corrected in [l5J. 
Refs. [ijjlljl provide analytical solutions valid for tensor 
modes that are super-horizon at decoupling, which give 
Ci spectra qualitatively valid for £ < 100. However this 
approximation was also used for / < 500 in Ref. |lj| i and 
so their result is qualitatively incorrect at I ^ 100. Their 
tensor polarization results also suffer qualitative prob- 
lems because the peak in the visibility was neglected. 
All previous analyses of sourced tensor modes have ne- 
glected neutrino compensation, giving results somewhat 
larger (but the result is still, in any case, somewhat un- 
certain) . 



VI. CONCLUSIONS 

Nearly scale invariant primordial magnetic fields can 
give a potentially observable CMB signal if they happen 
to be ^ lO'^^G. The observational B-mode signature 
comes from tensor modes, giving a non-Gaussian spec- 
trum otherwise essentially identical to that expected from 
primordial gravitational waves, and vector modes giving 
power on small scales. 

Any possible future detection of primordial gravita- 
tional waves from large scale B-modes should be carefully 
distinguished from that produced by magnetic fields. 
Any primordial signal is expected to be Gaussian, so 
Gaussianity tests can be used to distinguish them (meth- 
ods for robustly isolating the B-mode component on sec- 
tions of the sky are given in Ref. [23|, and can be used 
to construct a set of cut-sky modes that should be Gaus- 
sian if they are due to inflation). Low frequency obser- 
vations may also be able to detect Faraday rotation [l^ , 
which would be a clear signal of magnetic fields. Small 
scale B-mode observations from magnetic fields will need 
to be carefully distinguished from the weak lensing sig- 



nal. Regular primordial vector modes, though theoret- 
ically unmotivated, can also in principle give a signif- 
icant small and large scale B-mode signal [23 ■ They 
may be distinguished by their sharper fall in power on 
very small scales due to the lack of sources. Topological 
defects [^ can be identified by the lack of non-Gaussian 
tensor mode power on large scales. Note that throughout 
we have been assuming idealized observations; in practice 
foregrounds and systematics may well pose very serious 
problems (see e.g. Ref. "SS^). 

Our analysis is significantly more detailed than pre- 
vious work, in that we have numerically solved the full 
set of linearized equations. There is a qualitatively im- 
portant mechanism of a neutrino anisotropic stress com- 
pensation on super-horizon scales that was previously ne- 
glected. Computing the full transfer functions is rather 
straightforward, and we encourage future workers in this 
area to at least compare their semi-analytic results with 
the numerical answer to ensure that important physical 
effects have not be accidentally overlooked. Our modified 
version of CAMB^ [23 for efficiently computing vector 
mode power spectra is publicly available, and may also 
be useful for computing anisotropics from other sources, 
for example topological defects or second order effects. 
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^ An inconsistency between their definition Eq. 2.17 and the equa- 
tion for the C'l, Eq. 5.1 
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APPENDIX A: MULTIPOLE EQUATIONS, HARMONIC EXPANSION AND Ci 



In this appendix we review in a streamlined fashion the multipole equations, solutions, and equations for the C, 
needed for numerical calculation |24L l2a . l27l | . The definitions used here are precisely those used in the CAMB 
numerical implementation. Equivalent results using the total angular momentum method are given in Ref. p33 |. 

The photon multipole evolution is governed by the geodesic equation and Thomson scattering, giving [23 



IA: + gO/A, + D'hA, - ^[^D{alA,_,) + ^lAaJll - —laa.aM 



-Hear ( Iai - ISio - -;;IvaiSn - — Caiaj'^ia ) (AI) 



where Iai is taken to be zero for I < and 



Sab — . ^ab ~r ^t-'ab 



(A2) 



is a source from the anisotropic stress and E-polarization. The equation for the density perturbation Dal is obtained 
by taking the spatial derivative of the above equation for ^ = 0. The corresponding evolution equations for the 
polarization multipole tensors are |27| 



£a, + g©^^! + (1 + 1)2 ^^^bA, - ^^--pYD(a,fA,_i> - J-^ CuABa, = -TleariSAi - Y^Caia2^;2) 



^.,+^eS.,+ii±M_i)D^S,,,-^i^,,,^,,_, + ^curl£^ = 0. 



(A3) 



r 



For numerical solution we expand the covariant equa- 
tions into scalar, vector and tensor harmonics. The re- 
sulting equations for the modes at each wavenumber can 
be studied easily and also integrated numerically. 



2. Harmonic expansion 

For numerical work we perform a harmonic expansion 
in terms of zero order eigenfunctions of the Laplacian 



1. Scalar, vector, tensor decomposition 

It is useful to do a decomposition into ?Tj-type tensors, 
scalar (to — 0), vector (ttt, = 1) and 2-tcnsor (m = 2) 
modes. They describe respectively density perturbations, 
vorticity and gravitational waves. In general a rank— ^ 
PSTF tensor Xai can be written as a sum of ?7i— type 
tensors 



Xa. = E ^- 



,(m) 



Each Xj^ can be written in terms of I 
of a transverse tensor 



(A4) 



TO, derivatives 



D'Q"„ - ^Q!4„ 



(A6) 



where 






is transverse on all its indices, 
= 0. So a scalar would be expanded 
in terms of Q^, vectors in terms of Ql^, etc. We usually 
suppress the labelling of the different harmonics with 
the same eigenvalue, but when a function depends only 
on the eigenvalue we write the argument explicitly, e.g. 

fik). 

For m > there are eigenfunctions with positive and 
negative parity, which we can write explicitly as Q™ 
when required. Since 

1.2 

D2 ( curl QaJ= curl (D^ Q ^^ J = — curl Qa^ ( A7) 
(A5) they are related by the curl operation. Using the result 



where D, 



J-^ai J-^a2 



.D„ 



'Ai — J^ai J^a2 • • • J^a; 

PSTF and transverse D^-E 



and Sa„, is first order, 
A _ a — 0. The 'scalar' 

component is X'"-*, the 'vector' component is Xa , etc. 
Since General Relativity gives no sources for Xa^ with 
TO > 2 usually only scalars, vectors and (2-)tensors are 
considered. At linear order they evolve independently. 



curlcurlQ:i'„^ = — Q:^^ 



(A8) 



we can 


choose to normalize the ± harmonics the 


same 


way 


so 


that 














curlQ™± = 


S^A^- 


(A9) 



A rank-£ PSTF tensor of either parity may be con- 

")?n: 



structed from QT^ as 



S 



l~r 



D{A,_Q 



A^) 



(AlO) 



and an Xj^ component of Xai may be expanded in 
terms of these tensors. They satisfy 



D'Q% = Y^Q 



m 



2 ^2\ 






TO fc 



curlQ"± - --n™T 



(All) 



where the fc dependence of the harmonic coefficients is 
suppressed. We also often suppress m and ± indices for 
clarity. The other multipoles are expanded in analogy 
with Iai ■ The heat flux for each fluid component is given 
by Qi — {pi +Pi)vi, where u,; is the velocity, and the total 
heat flux is given by g = ^,- g^. We write the baryon ve- 



Jm) 



(m) 



locity simply as u, and define a constant Bq"^' = 11'^' /p^ 
to quantify the magnetic field anisotropic stress source. 
In the frame in which ^a = gradients are purely scalar 

(s'xf^ = 0. 



3. Harmonic multipole equations 



where I > m. 

Dimensionless harmonic coefficients are defined by 



Expanded into harmonics, the photon multipole equa- 
tions IjAip become 



"ab — Z_^ c'^ '^ab ^ ab ~ /^ c2 -" ^a 



k,± 



S 



k,± 



S^ 



b 



^W^^^Migjni Ei:^=J2LE( 



k,± 



k,± 



C2 " ^ab 



Jm) 



5:nM±g™^ /r = P.EA^ -A, 



7-(™)±/-)'"± 



fc,± 



k,± 



f^a = E ^^^Ql^ ^i"^ = E l^^-^^O"^ 



fc,± 



fc,± 



p,X)('")=Ee('5^)*"'^Qa 



rnib 



fc,± 



S*' 



(A12) 



II 



k 



2/ + 1 



(Z + l)2 



h+i — lli-i 



StleOT h - SiqIq - -SliV 



kaSi2 - 4:h'6u 



\kA5ii (A13) 



/ + 1 

2 
15' 

IS"""'" '■""" '3 
where I > m, Iq = Sp-y/p-y, Ii = for / < to, and to, su- 
perscripts are implicit. The scalar source is h' = {SS/S)' . 
The equation for the neutrino multipoles (after neutrino 
decoupling) is the same but without the Thomson scat- 
tering terms (for massive neutrinos see Ref. |39j). The 
polarization multipole equations l|A8|l become 



£™±' + k 



Bl''^' + k 



il + 3){l^ 


1)(/+1)2-to2 


{i + ir 

{l + 3){l- 


(2^ + 1) 

1)(? + 1)2-to2 



(i + ir 



(21 + 1) 



'-'l+l 



'-'l+l 



^ K?m± 

2Z-H1 '-1 



2m 



l{l + l 

2 TO 



/(/ + 1 



-K^ 



cm+ 
I i 



-Sn^ariSr^ - ^C^^n) 



(AM) 



4. Integral solutions 

Solutions to the Boltzmann hierarchies can be found in terms of line of sight integrals. The /™ hierarchy has 
homogeneous solutions (i.e. solutions to Eg. lAlSl with RHS set to zero) given by derivatives of 



*r(fc^) 



^! ii{kri) 
{I -my. {krj)"^ 



(A15) 



where ji{x) is a spherical Bessel function. These can be used to construct Green's function solutions to the full 
equations. For the polarization the result is less obvious, though solutions can easily be verified once found. For 
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vector modes (to = 1) the solutions are [32 



Mm) = 4 / drje '' 



Sne(JTV'i'l{x) + [ ka + Suear 



C\ d^lix) 



eHvo) 

B^im) 



i{i-i) r° 
T+rJ 

1(1-1) r^" 



drjSneCFTe ^ 
drjSneCfTe' 



l+l J ■ X 

where x = ^(^/o ~ v)- For tensors [ra = 2) the solutions are |27| 



1 djijx) ^ jijx) 
X dx X^ 



4/ dx 

± 



IiiVo) = 4 / drje-- 
1(1-1) 



no 



E^Vo) 
BHvo) = -2 



ka + SupCTT — 
4 



^?(x) 



(? + l)G + 2) 



drjSneCTe 



d^jiix) I 4 dj;(x) 



i-^b'(x) 



drjSnpdTe ^ 



(; + !)(/ + 2) 

Here t is the optical depth from 77 to 770, r' = —StieCFT 



djiix) , 2 
— ^ — + -J; X 
. c?x X 



(A16) 
(A17) 
(A18) 

(A19) 
(A20) 
(A21) 



5. Power spectra 

Using the harmonic expansion of Iai in Eq. 1181 the contribution to the Ci from type-?n tensors becomes 



a 



:TT 



(m) _ TT (2^! 



4(-2)'a!)\t7 



E (itkity)Q%uQy- 



(A22) 



The multipoles // can be related to some primordial variable Xa^ — X]fc {^^Q^a^ + ^ Q^a ) ^^^ ^ transfer function 
Tj"^ defined by /; = T(^ X . Statistical isotropy and orthogonality of the harmonics implies that 



{X^Xt,)^fx{k)5uu' 



(A23) 



where ^j, Skk'Yk = Yfe' and fx{k) is some function of the eigenvalue k. The normalization of the Q™ is given by 



k 



l—m 



dVQXQ-^' ^ Idvi 4Y"^ Q7^D^'-Q7, - i-2y-{i + my.il- my. ^ 



(20! 



(A24) 



where we have integrated by parts repeatedly, then repeatedly applied the identity for the divergence IjAlip . The 
normalization is A^ = J dVQ^"^QA^- By statistical isotropy Ci = (l/V) J dVCi and hence 



a 



TT{m) _TT [1 + my.{l - m)! Y^ N 



4 {-2r{l\)^ ^ V 



Y.i7\Ti''ik)\'fx{k) 



(A25) 



We choose to define a power spectrum Px{k) so that the real space isotropic variance is given by 

{\Xa^X^-\) = E ^/^W ^ j d\nk Px{k) (A26) 



A;,± 



so the CMB power spectrum becomes 



q"'"'=i "''g;i;);"'" /'"°^ftwi^."wi' 



(A27) 



Note that we have not had to choose a specific representation of (5a„ or ^j,. 
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The polarization Ci are obtained similarly ^7| and in general we have 



C, 



JK{m) 



G + l)(/ + 2) 



1(1-1) 



^1^ (20! (JmK^^) 



{-in 



H]\\^ 



p'. 



{l + l){l + 2)y''^ {I + m)\{l - m)\ 



1(1-1) 



2"(n)2 



I dink Px(k)J^Ki 



X r^X 



(A28) 
(A29) 



where JK is TT (p = 0), EE or BB (p = 2) or TE (p = 
1). Our conventions for the polarization are consistent 
with CMBFAST ^] and CAMB ^]. We have assumed 
a parity symmetric ensemble, so Cf^ 
For tensors we use Ht where hij = 



= cF^ = 0. 



Efc.±2i?Tg?, and 
hij is the transverse traceless part of the metric tensor. 
This introduces an additional factor of 1/4 into the result 
for the Ci in terms of Ph 



andT"^. 



The numerical factors in the hierarchy and Ci equa- 
tions depend on the choice of normalization for the £ and 
k expansions. Neither gm,) nor Qa, are normalized, so 
there are compensating numerical factors in the expres- 
sion for the Ci . If desired one can do normalized expan- 
sions, corresponding to an i!-dependent re-scaling of the 
Ii and other harmonic coefficients, givi ng expressions in 
more manifest agreement with Ref. |32l |. 
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